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I.  INTRODUCTION 

The  estimation  of  the  characteristics  of  the  ground 
motions  induced  by  near-surface  and  atmospheric  explosions  is 
a  complex  problem  which  has  received  extensive  study  for  more 
than  two  decades.  A  complete  solution  to  this  problem  would 
incorporate  an  explicit  consideration  of  the  highly  nonlinear 
response  of  the  medium  in  the  immediate  vicinity  of  ground 
zero,  including  the  vaporization  and  crushing  of  the  earth 
material,  leading  to  the  formation  of  a  surface  crater,  and 
the  subsequent  radiation  of  seismic  energy  away  from  the  cen¬ 
tral  source  region.  While  such  a  complete  solution  has  not 
yet  been  successfully  demonstrated,  it  has  recently  been  shown 
(Murphy,  1978;  Auld  and  Murphy,  1979;  Murphy  and  Bennett, 

1980)  that,  despite  these  complications,  a  significant  portion 
of  the  ground  motion  induced  by  surface  explosions  can  be  ex¬ 
plained  using  some  simple  linear,  elastic  approximations. 

In  particular,  it  has  been  concluded  that  the  late-time,  low- 
frequency  ground  motions  observed  from  such  explosions  are 
predominantly  elastic  Rayleigh  waves  produced  by  the  airblast 
loading  acting  on  the  surface  outside  the  region  of  strong 
nonlinear  response.  The  significance  of  this  finding  with 
regard  to  prediction  lies  in  the  fact  that  the  theory  govern¬ 
ing  the  generation  and  propagation  of  elastic  Rayleigh  waves 
is  well  developed  and  thus  can  be  used  to  help  define  scaling 
laws  for  the  low  frequency  components  of  the  ground  motion 
environments  produced  by  near-surface  explosions. 

The  research  summarized  in  this  report  has  focused  on 
the  extension  of  these  results  to  include  consideration  of 
atmospheric  explosions  detonated  over  a  wide  range  of  height 
of  burst  (HOB).  It  is  well  known  (  e.g.s  Brode ,  1968)  that 
the  airblast  associated  with  HOB  explosions  differs  in  many 
respects  from  that  associated  with  surface  explosions  of  the 
same  yield.  However,  up  until  now,  no  systematic  study  has 
been  made  of  how  this  variation  of  the  airblast  loading 


function  with  HOB  can  be  expected  to  affect  the  airblast- 
induced,  low  frequency  components  of  the  ground  motion.  An 
analytical  and  numerical  procedure  for  estimating  the  low 
frequency  ground  motions  produced  by  HOB  explosions  is  de¬ 
scribed  in  Section  II  where  the  results  of  a  sample  parame¬ 
tric  application  to  a  model  approximating  the  subsurface  geo¬ 
logy  at  Yucca  Flat  on  the  Nevada  Test  Site  is  also  presented. 
In  Section  III,  the  model  is  applied  to  the  prediction  of  the 
low  frequency  ground  motions  to  be  expected  from  the  Pre- 
Direct  Course  high  explosive  experiment  scheduled  to  be  con¬ 
ducted  at  the  White  Sands  Test  Site  in  the  fall  of  1982. 

This  is  followed  in  Section  IV  by  a  summary,  together  with  a 
statement  of  conclusions  regarding  the  current  state  of  knowl 
edge  with  respect  to  these  issues. 
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II.  RAYLEIGH  WAVES  FROM  NEAR-SURFACE  AND 
ATMOSPHERIC  EXPLOSIONS 

2.1  MODEL  DEFINITION 

Research  conducted  over  the  past  several  years  has 
provided  strong  evidence  that  the  low  frequency  components  of 
the  ground  motions  produced  by  near-surface  explosions  are 
predominantly  elastic  Rayleigh  waves  driven  by  the  airblast 
load  acting  on  the  surface  outside  the  region  where  strong 
nonlinear  effects  occur  (Murphy  and  Bennett,  1980;  Murphy, 
1981).  Thus,  the  boundary  value  problem  of  interest  for  pre¬ 
dicting  these  low  frequency  ground  motion  components  focuses 
on  estimation  of  the  elastic  response  of  a  multilayered  model 
of  the  subsurface  site  geology  to  a  propagating  normal  load 
acting  on  the  surface  (  i.e.,  an  approximation  to  an  airblast) 
The  geometry  of  the  problem  is  illustrated  in  Figure  1  which 
shows  a  vertical  section  consisting  of  an  axisymmetric  se¬ 
quence  of  constant  thickness  (d^)  layers  overlying  a  semi- 
infinite  halfspace  being  acted  on  by  a  concentrated  surface 
source.  As  a  first  approximation,  the  layers  are  taken  to  be 
perfectly  elastic  and  are  characterized  by  their  Lame  con¬ 
stants  (A,y)  and  density  (p).  Assuming  an  axisymmetric  air- 
blast  load,  the  Fourier  transform  of  the  Rayleigh  wave  com¬ 
ponent  of  the  displacement  for  a  particular  mode  at  observa¬ 
tion  point  r,  W(u,r),  can  be  written  as  a  spatial  integration 
over  the  Fourier  transform  of  the  overpressure,  p(r,w) ,  in 
the  form  (Murphy  and  Bennett,  1980) : 


W(w,r)  =  -iTrAR(oj)  Hq  (kr) 


p(r',w)  Jn  (kr ' )  r'  dr' 


+  J„(kr)  I  p(r',w)  h'  '  (kr')  r'  dr 


(1) 


where  AR(co)  is  the  Rayleigh  wave  site  response  function,  Jn , 

(2)  u 
Hq  denote  the  Bessel  and  Hankel  functions  respectively,  k 

is  the  wavenumber  of  the  Rayleigh  wave  component  with  fre¬ 
quency  gj  and  r^  denotes  the  exclusion  radius  which  defines 
the  region  within  which  strong  nonlinear  effects  occur.  For 
surface  bursts,  it  has  been  empirically  determined  that  Tq 
coincides  approximately  with  the  crater  radius.  However,  for 
atmospheric  explosions  with  significant  HOB  (  i.e.  ,  greater 
than  Tq) ,  it  seems  likely  that  this  heuristic  approximation 
can  be  avoided  and  rg  set  equal  to  zero.  Thus,  comparisons 
of  predicted  and  observed  surface  waves  for  explosions  at 
significantly  different  HOB  can  provide  a  definitive  test  of 
the  adopted  source  approximation  for  near-surface  bursts. 

2.2  AIRBLAST  LOADING  FROM  ATMOSPHERIC  EXPLOSIONS 

The  analytic  surface  wave  simulations  conducted  to 
date  have  dealt  exclusively  with  surface  explosions,  and  the 
extension  to  atmospheric  explosions  with  arbitrary  HOB  re¬ 
quires  some  modifications  in  the  prediction  procedure.  In 
particular,  although  equation  (1)  is  still  valid,  a  modified 
expression  is  required  for  the  airblast  loading,  p(r,t),  and 
its  corresponding  Fourier  transform,  p(r,co).  For  this  pur¬ 
pose,  we  have  adopted  Brode ' s  recently  revised  approximation 
(Speicher  and  Brode,  1981)  which  expresses  the  overpressure 
for  a  given  yield  (W)  and  HOB  as  a  function  of  range  and  time 
through  products  of  functions  involving  fairly  elaborate 
polynomial  fits  to  the  characteristics  of  airblasts  observed 
from  surface  and  atmospheric  nuclear  explosions.  The  princi¬ 
pal  features  of  the  overpressure  time  histories  predicted  by 
this  model  are  illustrated  schematically  in  Figure  2.  Note 
that  the  model  predicts  only  the  positive  phase  of  the  over¬ 
pressure  beginning  at  the  shock  arrival  time  t  and  having 
duration  Dp.  Also,  in  contrast  to  the  overpressure  observed 
from  surface  bursts  which  monotonically  decay  from  a  peak 


value  at  the  shock  arrival  time  t  «  t  ,  the  overpressure  from 
HOB  explosions  can  show  significant  secondary  peaks,  such  as 
that  indicated  in  Figure  2 ,  as  a  result  of  the  Mach  reflection 
of  the  incident  shock  wave  from  the  free  surface.  Thus,  in 
general,  finer  time  resolution  is  required  to  adequately  de¬ 
fine  the  characteristics  of  the  airblast  loading  induced  by 
HOB  explosions . 

Some  of  the  prominent  characteristics  of  the  airblast 
loading  as  a  function  of  HOB  at  W  =  1  kt  are  illustrated  in 
Figures  3-6.  Figure  3  compares  the  peak  incident  pressure, 
p(t_) ,  as  a  function  of  range  for  HOB  of  0,  100  and  500  m. 

a. 

As  might  be  expected,  at  ranges  less  than  the  HOB,  the  pre¬ 
dicted  peak  pressures  for  atmospheric  explosions  are  consid¬ 
erably  lower  than  those  predicted  for  a  corresponding  surface 
burst.  However,  the  low  frequency  components  of  the  induced 
ground  motions  are  relatively  insensitive  to  variations  in 
peak  pressure  and  it  is  necessary  to  consider  how  the  other 
characteristics  of  the  overpressure  time  histories  vary  with 
HOB.  Figure  4  compares  shock  arrival  time,  t  ,  and  duration 
of  positive  phase,  D  ,  as  a  function  of  range  for  a  surface 
burst  and  an  atmospheric  explosion  with  HOB  =  500  m.  Again, 
for  ranges  less  than  the  HOB,  the  Dp  values  for  the  two  cases 
differ  significantly,  indicating  changes  in  the  time  history 
of  the  airblast  loading  in  the  vicinity  of  ground  zero.  This 
is  illustrated  more  specifically  in  Figure  5  which  compares 
the  overpressure  time  histories  as  a  function  of  HOB  at  a 
range  of  100  m.  It  can  be  seen  that  although  the  peak  pres¬ 
sure  decreases  with  increasing  HOB,  as  predicted  in  Figure  3, 
this  effect  on  the  positive  phase  impulse  is  largely  offset 
by  the  increase  in  Dp  with  HOB.  The  corresponding  overpres¬ 
sure  time  histories  as  a  function  of  HOB  at  a  range  of  1000  m 
are  compared  in  Figure  6.  This  figure  illustrates  that  at 
ranges  greater  than  the  largest  HOB,  the  airblast  loading  is 
relatively  insensitive  to  HOB. 


The  results  shown  in  Figures  3-6  indicate  that  the  air 
blast  loading  shows  a  complicated  dependence  on  HOB  which  can 
be  expected,  in  general,  to  have  corresponding,  frequency- 
dependent  effects  on  the  induced  surface  waves.  However,  as 
has  been  noted  previously  (Murphy  and  Bennett,  1980),  some 
simplifications  are  possible  in  the  low  frequency  limit. 

That  is,  as  w  — »  0 ,  the  amplitude  of  the  airblast-induced  Ray¬ 
leigh  waves  are  directly  proportional  to  the  total  impulse 
delivered  by  the  airblast  over  the  loaded  area.  Figure  7 
shows  the  total  impulse  as  a  function  of  HOB  delivered  by  a 
1  kt  explosion  inside  a  radius  of  2  km  from  ground  zero 
the  radius  at  which  the  peak  overpressure  drops  to  a  level  of 
about  1  psi) .  It  can  be  seen  that  this  total  impulse  varies 
by  less  than  a  factor  of  1.5  for  HOB  extending  from  0  to  500 
m.  Thus,  despite  the  orders  of  magnitude  variations  in  peak 
pressure  shown  in  Figure  3,  the  total  impulse  is  relatively 
insensitive  to  HOB,  which  suggests  that  the  low  frequency 
components  of  the  induced  ground  motions  will  not  be  very 
strongly  dependent  on  HOB.  This  dependence  will  be  explored 
more  explicitly  in  the  following  Sections  2.3  and  III. 

For  applications  employing  equation  (1);  the  Fourier 
transform  of  the  overpressure,  p(r,a>),  is  required  and  it  is 
not  feasible  to  compute  this  analytically  from  the  p(r,t) 
given  by  the  complex  polynomial  approximations  of  Speicher 
and  Brode  (1981).  Thus,  a  numerical  Fourier  transform  proce¬ 
dure  must  be  employed.  A  problem  arises,  however,  in  that 
the  overpressure  time  histories  are  often  only  a  fraction  of 
a  second  in  total  duration  (cf.  Figure  5).  Thus,  the  lowest 
nonzero  frequency  provided  by  a  numerical  Fourier  transform 
routine  will  be  on  the  order  of  10  Hz,  while  the  interest  for 
surface  wave  simulation  purposes  is  almost  always  in  compo¬ 
nents  of  the  motion  with  frequencies  less  than  a  few  Hertz. 
Thus,  a  procedure  is  required  for  interpolating  the  numerical 
Fourier  transform  values  in  the  frequency  range  of  interest. 


This  problem  has  been  extensively  studied  (  e.g Papoulis, 
1962)  and  the  solution  for  the  present  application  may  be 
summarized  as  follows.  Suppose  the  overpressure  at  a  given 
range,  p(x),  is  defined  such  that  p(x)  =  0  for  t  <  0  and  x  > 
2T  with  x  =  t-t  .  Then,  it  follows  that  its  Fourier  trans- 

a 

form  can  be  expressed  in  the  form 

J2T  -iojx 

p(x)  e  dx  (2) 

0 

Then,  by  definition 

.  nr 

r  2T  1  T  T 

Pn  E  P  p{-)  =  J  p(t)  e  dx  (3) 

0 


On  the  other  hand,  if  p(x)  is  expanded  in  a  Fourier  series  in 
the  interval  0  <  x  <  2T,  one  obtains 


P (T )  =  £ 


an  e 


TT 

in  j  t 


n=  -° 


where 


•  j_  f  2T 

an  2T  j  P 


-in  ip  x 


(x)  e 


dx 


It  follows  from  (3)  and  (5)  that 


n 


n 


2T 


(4) 


(5) 


(6) 


and  a  substitution  of  (6)  and  (4)  into  (2)  gives 


T 


dx 


sin  2  (ooT-nTr ) 
(oiT  -  nir) 


.  2sin  (oiT-niT) 
1  (ojT  -  nTT ) 


and  thus  we  can  write 


p(co) 


-nir . 
P(xp-) 


sin  (coT  -  nir)  _i  (coT  “  n1T) 
(coT  -  rnr)  e 


Finally,  using  the  identity 


inir 

e  sin(wT  -  nir)  =  sin  coT 


we  obtain 


p(co) 


-ia)T 

e  sin  coT 

n=  -oo 


,  nn . 
P(-f-) 

(coT  -  niT) 


(1 


(1 


Setting  p  =  pR  +  i  p_,  the  interpolation  algorithms  for  the 
real  and  imaginary  parts  of  the  Fourier  transform  of  the 
pressure  function  take  the  form 


PRU) 


=  sin  cuT 


n=  -o 


=  sin  coT 


Pr(tjt)  cos  coT  +  pT(^?-)  sin  coT 


I '  T 


(wT  -  nxT ) 


(1 


/  nTT  x  _  /  n TT  x  .  _ 

PI(-^-)  cos  0)T  -  Pr(-^t)  Sin  ojT 


d,  (co ) 


which  expresses  the  Fourier  transform  of  the  overpressure  at 
an  arbitrary  frequency  co  in  terms  of  the  evenly  spaced  values 
obtained  from  the  numerical  Fourier  transform.  In  actual 
applications,  the  pressure  at  a  given  range  is  nonzero  over 
the  time  interval  extending  from  ta  to  tfl  +  Dp.  Thus,  we  set 
T  =  D  /2  in  (12)  and  then  the  required  values  of  the  real  and 

tr  ^  ^ 

imaginary  parts  of  the  pressure  (pR,  Pj) >  referenced  to  real 
time  t,  are  given  by 

pR(co)  =  pR(co)  cos  cot  +  Pj(co)  sin  cota 

(13) 

Pj(a)  =  p-j-(co)  cos  cota  -  pR(co)  sin  cota 

2.3  SIMULATED  RAYLEIGH  WAVES  FROM  EXPLOSIONS  WITH 

DIFFERENT  HOB 

In  order  to  illustrate  the  simulation  of  Rayleigh 
waves  as  a  function  of  HOB,  we  will  consider  a  preliminary 
application  to  a  subsurface  geologic  model  characteristic  of 
the  Yucca  Flat  testing  area  on  the  Nevada  Test  Site.  Esti¬ 
mates  of  the  relevant  subsurface  geophysical  parameters  (com- 
pressional  and  shear  wave  velocities  a  and  B  and  density  p) 
for  this  site  are  illustrated  graphically  in  Figure  8  (Hays 
and  Murphy,  1971)  .  It  can  be  seen  that  the  area  is  charac¬ 
terized  by  thick  layers  of  alluvium  and  tuff  overlying  the 
Paleozoic  basement  rocks.  The  fundamental  mode  Rayleigh  wave 
dispersion  curves  computed  using  this  model  are  shown  in 
Figure  9,  where  it  can  be  seen  that  there  is  a  well-developed 
group  velocity  minimum  in  the  0.5  to  1.0  Hz  frequency  band. 

As  has  been  noted  previously  (Murphy  and  Bennett,  1980),  this 
frequency  of  minimum  group  velocity  coincides  with  a  corner 
frequency  in  the  Rayleigh  wave  site  response  function  and  is 
given  approximately  by  8/2H  where  H  denotes  the  depth  to  the 
basement  and  B  denotes  the  interval-averaged  shear  wave  veloc¬ 
ity  to  that  depth.  This  is  illustrated  in  Figure  10  which 
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Figure  8 


Alluvium 

a  =  1600  m/sec 
B  =  1000  m/sec 
p  =  1 . 9  gm/cm^ 


Tuff 

a  =  2650  m/sec 
B  =  1550  m/sec 
p  =  1.9  gm/cm^ 
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Paleozoics 
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B  =  3000  m/sec 

3 

p  =  2.8  gm/cm 


* 

Subsurface  geologic  model  for  Yucca  Flat  simula¬ 
tions  . 
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shows  the  corresponding  Rayleigh  wave  site  response  function 
with  a  pronounced  corner  frequency  at  about  0.7  Hz.  For  ex¬ 
plosions  of  sufficiently  large  yield,  this  corner  frequency 
will  define  the  predominant  frequency  of  the  explosion-induced 
Rayleigh  waves  (Murphy  and  Bennett,  1980). 

The  fundamental  mode  Rayleigh  waves  induced  in  this 
geologic  model  by  1  kt  explosions  with  HOB  ranging  from  0  to 
500  m  have  been  computed  using  equation  (1)  together  with  the 
airblast  approximation  of  Speicher  and  Brode  (1981)  .  The  re¬ 
sulting  Rayleigh  wave  displacement  spectra  at  a  range  of 
2500  m  from  1  kt  explosions  at  HOB  of  100  and  500  m  are  com¬ 
pared  in  Figure  11.  It  can  be  seen  that  the  high  frequency 
content  of  the  signal  decreases  significantly  with  increasing 
HOB  over  this  range.  This  presumably  reflects  the  relatively 
low  frequency  character  of  the  incident  overpressure  near 
ground  zero  for  the  higher  altitude  explosion  (cf .  Figure  3) . 
However,  these  spectral  differences  are  most  notable  in  a 
frequency  band  within  which  the  spectral  amplitudes  are  an 
order  of  magnitude  down  from  the  peak  level  and,  thus,  con¬ 
tribute  little  of  the  time  domain  signal.  In  the  band  be¬ 
tween  0.6  and  1.0  Hz,  in  which  the  spectral  maxima  occur,  the 
differences  between  the  two  spectra  are  much  less  pronounced. 
Moreover,  it  has  been  found  that  for  HOB  in  the  range  0  to 
100  m,  the  Rayleigh  wave  displacement  spectra  are  essentially 
identical  to  that  obtained  for  an  HOB  of  100  m.  Thus,  the 
spectral  results  suggest  that  the  airblast-induced  Rayleigh 
waves  at  this  site  should  be  relatively  insensitive  to  HOB. 
This  inference  is  confirmed  by  Figures  12-14  which  show  the 
dependence  of  the  induced  Rayleigh  waves  on  HOB  at  ranges  of 
500,  1000  and  2500  m,  respectively.  Comparing  Figures  12-14 
for  a  particular  value  of  HOB  it  can  be  seen  that,  as  ex¬ 
pected,  the  dominant  frequency  of  the  Rayleigh  waves  at  all 
three  ranges  fall  in  the  0.7  to  1.0  Hz  band  and  the  predomi¬ 
nant  effect  of  increasing  range  on  the  signal  is  the  gradual 
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Figure  12.  Vertical  component  Rayleigh  wave  displacements 
as  a  function  of  HOB  at  a  range  of  500  m,  Yucca 
Flat,  W  =  1  kt. 
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Figure  13.  Vertical  Component  Rayleigh  wave  displacement 
as  a  function  of  HOB  at  a  range  of  1000  m,  Yu< 
Flat,  W  =  1  kt. 


t,  seconds 


Figure  14.  Vertical  component  Rayleigh  wave  displacements 

as  a  function  of  HOB  at  a  range  of  2500  m,  Yucca 
Flat,  W  =  1  kt. 


dispersion  of  energy  over  time  due  to  the  frequency  dependent 
propagation  velocities  of  the  various  components  of  the  signal 
Furthermore,  it  can  be  seen  that  at  a  fixed  range,  the  simu¬ 
lated  Rayleigh  waves  show  little  variation  with  HOB  over  the 
range  from  0  to  500  m.  Specifically,  the  frequency  content 
remains  essentially  constant  and  the  amplitude  level  varies 
by  less  than  a  factor  of  two.  Thus,  as  was  inferred  from  the 
integrated  impulse  curve  shown  in  Figure  7,  the  low  frequency 
airblast-induced  Rayleigh  waves  are  relatively  insensitive  to 
HOB,  at  least  for  this  combination  of  site  geology  and  yield. 


III.  APPLICATION  TO  THE  PRE-DIRECT  COURSE  EXPERIMENT 


The  Pre-Direct  Course  high  explosive  (HE)  experiment 
is  scheduled  to  be  detonated  at  the  White  Sands  Test  Site  in 
the  fall  of  1982,  near  the  location  of  the  previous  Mill  Race 
experiment.  As  currently  planned,  it  will  consist  of  24  tons 
of  HE  (ANFO)  detonated  at  an  HOB  of  16.7  m.  Preliminary 
determinations  of  subsurface  physical  properties  at  the  site 
have  been  made  on  the  basis  of  a  refraction  survey  (Personal 
Communication,  Bob  Reinke,  AFWL ,  1982)  and  are  shown  schemat¬ 
ically  in  Figure  15  together  with  the  corresponding  predicted 
fundamental  mode  Rayleigh  wave  dispersion  curves.  It  should 
be  noted  that  the  shear  wave  velocities,  which  control  the 
Rayleigh  wave  characteristics ,  have  not  been  measured  at  this 
site  and,  consequently,  have  been  estimated  on  the  basis  of 
results  obtained  from  other  sites  in  this  general  area.  It 
can  be  seen  that  the  refraction  survey  indicates  the  presence 
of  a  thick  (75m)  layer  of  dry  alluvium  overlying  saturated 
alluvium  which  apparently  extends  to  a  depth  of  more  than 
150  m  (  i.e. ,  the  approximate  limit  of  resolution  of  the  500  m 
long  refraction  line) .  Note  that  because  of  the  small  in¬ 
ferred  variation  of  shear  wave  velocity  (6)  with  depth,  very 
strong  dispersion  of  the  surface  wave  energy  is  not  expected 
at  this  site.  However,  the  model  does  predict  a  distinct 
group  velocity  minimum  at  a  frequency  of  about  2 . 5  Hz  (  i.e., 
Bj/2H)  .  As  expected,  this  corresponds  to  a  pronounced  corner 
frequency  on  the  corresponding  Rayleigh  wave  site  response 
function,  Ad(oj),  shown  in  Figure  16.  Thus,  for  this  site 
model,  the  predominant  energy  in  the  surface  wave  components 
of  the  ground  motion  is  expected  to  be  confined  to  frequencies 
above  2.5  Hz . 

The  low  frequency  ground  motions  to  be  expected  from 
this  test  have  been  estimated  using  an  airblast  loading  func¬ 
tion  computed  from  the  Speicher  and  Brode  (1981)  analytic 


approximation  for  nuclear  explosions  with  an  effective  yield 
of  48  tons.  Available  evidence  suggests  that  a  nuclear  explo¬ 
sion  of  this  yield  will  produce  an  airblast  loading  essen¬ 
tially  equivalent  to  that  expected  from  the  detonation  of  24 
tons  of  ANFO  at  that  same  height  (Personal  Communication, 
George  Ullrich,  DNA,  1982).  Figure  17  shows  a  comparison  of 
the  predicted  peak  overpressures  as  a  function  of  range  ob¬ 
tained  for  Pre-Direct  Course  and  a  surface  burst  of  the  same 
yield.  The  surface  burst  prediction  is  truncated  at  a  range 
of  6  m,  consistent  with  the  estimated  exclusion  radius  (  i.e., 
Tq)  for  that  yield.  It  can  be  seen  that  the  predicted  peak 
pressures  are  quite  different  at  ranges  of  less  than  15  m  and 
essentially  identical  at  larger  ranges. 

The  airblast  induced  Rayleigh  waves  to  be  expected 
from  Pre-Direct  Course,  as  well  as  from  an  expected  subsequent 
companion  surface  burst  experiment  of  the  same  yield,  have 
been  estimated  using  the  analytical  models  and  numerical  pro¬ 
cedures  described  in  Section  II.  Figure  18  shows  the  pre¬ 
dicted  vertical  component  Rayleigh  wave  displacements  for 
Pre-Direct  Course  over  the  horizontal  distance  range  extend¬ 
ing  from  ground  zero  to  200  m.  It  can  be  seen  that  the  pre¬ 
dicted  waveforms  are  quite  simple  and  pulse-like,  reflecting 
the  simplicity  of  the  adopted  subsurface  geological  model 
{i.e.,  the  shear  wave  velocity  distribution  shown  in  Figure 
15  is  nearly  independent  of  depth) .  The  most  notable  feature 
of  these  waveforms  is  the  evidence  of  multiple  arrivals  at 
ranges  greater  than  100  m.  It  was  initially  conjectured  that 
this  multiple  arrival  phenomenon  might  be  associated  with  the 
secondary  Mach  reflection  arrival  in  the  airblast  loading 
predicted  for  the  HOB  explosion.  However,  it  was  subsequently 
noted  that  the  delay  time  between  the  two  airblast  arrivals 
is  much  too  short  to  be  associated  with  the  observed  delay 
between  the  multiple  Rayleigh  wave  arrivals.  In  fact,  the 
results  of  numerical  modeling  experiments  suggest  this 


33 


0  0.5  1.0  1-5 

Time,  seconds 


Figure  18.  Predicted  vertical  component  Rayleigh  wave  dis¬ 
placement  as  a  function  of  range,  Pre-Direct 
Course . 


characteristic  of  the  predicted  Rayleigh  waves  is  related  to 
the  characteristics  of  the  site.  That  is,  the  Rayleigh  wave 
coupling  for  an  explosion  of  this  yield  for  this  site  model 
is  partitioned  into  two  distinct  phases:  an  initial  phase 
dominated  by  the  loading  within  about  20  m  of  ground  zero  and 
a  secondary  phase  dominated  by  the  loading  in  the  range  from 
about  50  to  100  m.  This  is  presumably  a  function  of  the  re¬ 
lationship  between  the  local  airblast  velocity  and  the  Ray¬ 
leigh  phase  velocity  distribution  characteristic  of  the  site. 
Thus,  the  predicted  effect  is  probably  strongly  dependent  on 
the  subsurface  shear  wave  velocity  distribution,  which  is  not 
very  well  constrained  by  currently  available  data. 

As  in  the  case  of  the  Yucca  Flat  nuclear  simulations 
described  in  Section  2.3,  it  has  been  found  that  the  predicted 
Rayleigh  displacements  for  an  equivalent  surface  burst  test 
are  quite  similar  to  those  obtained  for  Pre-Direct  Course. 

This  is  illustrated  in  Figure  19  which  shows  comparisons  of 
the  predicted  Rayleigh  waves  for  the  two  explosions  at  ranges 
of  25  and  200  m.  It  can  be  seen  that  the  predicted  peak  dis¬ 
placements  and  waveforms  agree  to  within  the  expected  measure¬ 
ment  uncertainty  and  indicate  no  pronounced  effect  of  HOB 
over  this  range. 

Thus,  the  theoretical  simulations  of  the  Rayleigh  waves 
to  be  expected  from  the  Pre-Direct  Course  HOB  experiment  and 
a  surface  burst  of  the  same  yield  suggest  that  there  should 
be  negligible  differences  in  the  low  frequency  ground  motions 
observed  from  these  two  tests.  It  follows  that  ground  motion 
data  recorded  from  such  a  pair  of  experiments  conducted  at 
this  site  would  provide  an  excellent  test  of  the  validity  of 
the  approximations  currently  being  used  to  estimate  this 
component  of  the  explosively-generated  ground  motions. 


IV.  SUMMARY  AND  CONCLUSIONS 


4 . 1  SUMMARY 

The  investigation  summarized  in  this  report  has  cen¬ 
tered  on  an  analysis  of  the  effects  of  variations  of  HOB  on 
the  low  frequency  ground  motions  induced  by  near-surface  and 
atmospheric  explosions.  In  particular,  the  present  effort 
has  encompassed  an  extension  of  the  existing  analytical  model 
for  surface  explosions  to  include  consideration  of  a  general 
airblast  loading  function  valid  for  arbitrary  HOB. 

The  mathematical  model  derived  to  estimate  the  low 
frequency  ground  motions  produced  by  HOB  explosions  was  de¬ 
scribed  in  Section  II  where  it  was  applied  to  the  parametric 
investigation  of  the  effects  of  HOB  on  a  prototype  1  kt 
nuclear  explosion  detonated  over  a  site  model  approximating 
the  subsurface  geology  at  Yucca  Flat  on  the  Nevada  Test  Site. 
The  generalized  airblast  loading  function  proposed  by  Speicher 
and  Brode  (1981)  was  adopted  for  modeling  purposes  and  a  de¬ 
tailed  numerical  procedure  was  derived  for  incorporating  this 
loading  function  into  the  analytical  Rayleigh  wave  prediction 
model.  The  results  of  the  parametric  investigation  of  the 
Yucca  Flat  site  model  were  analyzed  in  both  the  frequency  and 
time  domains  and  it  was  found  that,  despite  the  fact  that  the 
characteristics  of  the  airblast  loading  changes  significantly 
with  increasing  HOB,  the  corresponding  induced  surface  waves 
are  relatively  insensitive  to  HOB.  That  is,  it  was  found 
that  the  low  frequency  coupling  of  interest  with  regard  to 
surface  wave  generation  is  not  very  sensitive  to  HOB,  at 
least  for  this  combination  of  site  geology  and  yield. 

The  estimation  of  the  characteristics  of  the  low  fre¬ 
quency  ground  motions  to  be  expected  from  the  Pre-Direct 
Course  HOB  experiment  was  addressed  in  Section  III  where 
predictions  of  near- field  ground  motions  were  presented  for 
this  event  as  well  as  for  a  hypothetical  surface  burst  test 


of  the  same  yield  at  that  site.  It  was  shown  that  the  surface 
waves  predicted  on  the  basis  of  the  currently  available  sub¬ 
surface  model  of  the  site  are  quite  simple  and  pulse-like  and 
show  little  dependence  on  HOB  over  the  range  considered. 

4.2  CONCLUSIONS 

The  analyses  summarized  above  support  the  following 
conclusions  regarding  the  effects  of  HOB  on  the  low  frequency 
ground  motions  induced  by  near-surface  and  atmospheric  explo¬ 
sions  . 

1.  Despite  the  fact  that  the  airblast  loading 
from  explosions  varies  significantly  with 
HOB,  the  total  integrated  impulse  delivered 
to  the  ground  by  the  propagating  airblast, 
which  defines  the  low  frequency  limit  of  the 
surface  wave  coupling  coefficient,  is  rela¬ 
tively  insensitive  to  HOB  for  scaled  HOB  in 
the  range  from  0  to  500  m/kt^^. 

2.  Simulated  Rayleigh  waves  from  1  kt  explosions 
over  a  site  model  approximating  the  subsurface 
geology  at  Yucca  Flat  on  the  Nevada  Test  Site 
were  found  to  be  essentially  independent  of 
HOB  for  HOB  ranging  from  0  to  500  m.  This  is 
consistent  with  Cooper's  (1972)  observation 
that  the  peak  displacement  data  recorded  from 
nuclear  airblasts  at  NTS  are  not  strongly 
dependent  on  HOB. 

3.  Predicted  surface  wave  ground  motions  for 
the  proposed  Pre-Direct  Course  HE  experiments 
and  a  hypothetical  surface  burst  test  of  the 
same  yield  also  show  little  dependence  on  HOB. 

For  this  reason,  it  is  concluded  that  such  a 
pair  of  experiments  would  provide  an  excellent 
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Naval  Civil  Engineering  Laboratory 
ATTN:  Code  LSI,  J.  Crawford 

Naval  Electronic  Systems  Command 
ATTN:  PME  117-21 

Naval  Facilities  Engineering  Command 
ATTN :  Code  04B 

Naval  Material  Command 
ATTN:  MAT  08T-22 

Naval  Postgraduate  School 
ATTN:  G.  Lindsay 
ATTN:  Code  1424  Library 

Naval  Research  Laboratory 
ATTN:  Code  2627 


Naval  Surface  Weapons  Center 
ATTN:  Code  F31 

Naval  Surface  Weapons  Center 

ATTN:  Tech  Library  &  Info  Svcs  Br 

Naval  War  Col  lege 

ATTN:  Code  E-ll,  Tech  Service 


43 


DEPARTMENT  OF  THE  NAVY  (Continued) 


Naval  Weapons  Evaluation  Facility 
ATTN:  Code  10 


Newport  Laboratory 

ATTN:  Code  EM,  J.  Kalinowski 


Ofc  of  the  Deputy  Chief  of  Naval  Ops 
ATTN:  OP  03EG 
ATTN:  NOP  981 


Office  of  Naval  Research 

ATTN:  Code  474,  N.  Perrone 


Strategic  Systems  Project  Office 
ATTN:  NSP-43 


DEPARTMENT  OF  THE  AIR  FORCE 


Air  Force 
ATTN: 


INT 


Air  Force  Geophysics  Laboratory 
ATTN:  LWH,  H.  Ossing 


Air  Force  Institute  of  Technology 
ATTN:  Library 


Air  Force  Office  of  Scientific  Rsch 
ATTN:  W.  Best 
ATTN:  J.  Allen 


Air  Force  Systems  Command 


Deputy  Chief  of  Staff 
Logistics  &  Engineering 
ATTN:  LEEE 


Deputy  Chief  of  Staff 
Research,  Development  &  Acq 
ATTN:  AFRDQI 


Foreign  Technology  Division 
ATTN:  NIIS  Library 


Rome  Air  Development  Center 
ATTN:  TSLD 


Strategic  Air  Command 

ATTN:  INT,  E.  Jacobsen 
ATTN:  NRI-STINFO  Library 


DEPARTMENT  OF  ENERGY 


Albuquerque  Operations  Office 
ATTN:  R.  Jones 
ATTN:  CTID 


Nevada  Operations  Office 

ATTN:  Doc  Con  for  Technical  Library 


OTHER  GOVERNMENT  AGENCIES 


Central  Intel  1 igence  Agency 
ATTN:  OSWR/NED 


Department  of  the  Interior 
Bureau  of  Mines 

ATTN:  Tech  Lib 


NATO 


NATO  School  (SHAPE) 

ATTN:  US  Documents  Officer 


DEPARTMENT  OF  ENERGY  CONTRACTORS 


University  of  California 
Lawrence  Livermore  National  Lab 
ATTN:  D.  Glenn 

ATTN:  Technical  Info  Dept  Library 


Los  Alamos  National  Laboratory 
A1TN:  Reports  Library 
ATTN:  J.  Hopkins 


ATTN:  DLW 

Oak  Ridge  National  Laboratory 

ATTN:  Civil  Def  Res  Proj 

Air 

Force  Weapons  Laboratory 

ATTN:  SUL 

ATTN:  Central  Rsch  Library 

ATTN:  NTES,  J.  Thomas 

Sandia  National  Lab 

ATTN:  NTE,  M.  Plamondon 

ATTN:  Tech  Lib,  3141 

ATTN:  NTES-C,  R.  Henny 

ATTN:  Org  7112,  A.  Chabai 

ATTN:  L.  Hill 

Air 

University  Library 

ATTN:  AUL-LSE 

Sandia  National  Labs,  Livermore 

ATTN:  Library  &  Security  Classification 

Ballistic  Missile  Office 

ATTN:  ENBF,  D.  Gage 

ATTN:  ENSN,  A.  Schenker 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ATTN:  EN 

Aerospace  Corp 

ATTN:  Technical  Information  Services 

-4 


Agbabian  Associates 

ATTN:  M.  Agbabian 


Applied  Research  Associates,  Inc 


ATTN 

ATTN 

ATTN 


N.  Higgins 
J.  Bratton 
H.  Auld 


Applied  Research  Associates,  Inc 
ATTN:  S.  Blouin 


Applied  Research  Associates,  Inc 
ATTN:  D.  Piepenburg 


Applied  Research  Associates,  Inc 
ATTN:  B.  Frank 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued 


Applied  Theory,  Inc 
2  cy  ATTN:  J.  Trulio 


Lockheed  Missiles  &  Space  Co,  Inc 

ATTN:  Technical  Information  Center 
ATTN:  T.  Geers 


AVCO  Systems  Division 


ATTN: 

Library  A830 

Lockheed  Missiles  &  Space  Co,  Inc 

ATTN 

J.  Weisner,  Dept  80-82 

BDM  Corp 

ATTN 

T.  Neighbors 

Martin  Marietta  Denver  Aerospace 

ATTN 

Corporate  Library 

ATTN 

D-6074,  G.  Freyer 

BDM  Corp 

McDonnell  Douglas  Corp 

ATTN 

F.  Leech 

ATTN 

R.  Halprin 

Boeing  Aerospace  Co 

Merritt  CASES,  Inc 

ATTN 

MS-85-20,  D.  Choate 

ATTN 

Library 

ATTN 

J.  Merritt 

Boeing  Co 

ATTN 

Aerospace  Library 

University  of  New  Mexico 

ATTN 

N.  Baum 

Cal ifornia 

Institute  of  Technology 

ATTN 

T.  Ahrens 

Pacific-Sierra  Research  Corp 

ATTN 

H.  Brode,  Chairman  SAGE 

Cal ifornia 

Research  &  Technology,  Inc 

ATTN 

S.  Schuster 

Patel  Enterprises,  Inc 

ATTN 

Library 

ATTN 

M.  Patel 

ATTN 

K.  Kreyenhagen 

ATTN 

M.  Rosenblatt 

Physics  International  Co 

ATTN 

J.  Thomsen 

California 

Research  &  Technology,  Inc 

ATTN 

F.  Sauer 

ATTN 

D.  Orphal 

ATTN 

E.  Moore 

ATTN 

Technical  Library 

Cal  span  Corp 

ATTN 

L.  Behrmann 

ATTN 

Library 

R&D  Associates 

University 

of  Denver 

ATTN 

R.  Port 

ATTN:  Sec  Officer  for  J.  Wisotski 

ATTN 

J.  Carpenter 

ATTN 

W.  Wright 

EG&G  Wash  Analytical  Svcs  Ctr,  Inc 

ATTN 

Technical  Information  Center 

ATTN 

Library 

ATTN 

J.  Lewis 

Electro-Mech  Systems,  Inc 

R&D  Associates 

ATTN 

R.  Shunk 

ATTN 

H.  Cooper 

Gard,  Inc 

S-CUBED 

ATTN 

G.  Neidhardt 

ATTN 

Library 

ATTN 

D.  Grine 

Horizons  Technology,  Inc 

ATTN 

T.  Riney 

ATTN:  R.  Kruger 

ATTN 

T.  Cherry 

4  cy  ATTN 

J.  Murphy 

I IT  Research  Institute 

4  cy  ATTN 

H.  Shah 

ATTN 

M.  Johnson 

4  cy  ATTN 

T.  Tzeng 

ATTN 

R.  Welch 

ATTN 

Documents  Library 

S-CUBED 

ATTN 

J.  Murphy 

Institute  for  Defense  Analyses 

ATTN:  Classified  Library 

Science  Applications,  Inc 

ATTN 

Technical  Library 

Kaman  AviDyne 

ATTN 

Library 

Science  Appl ications,  Inc 

ATTN 

N.  Hobbs 

ATTN 

D.  Maxwell 

ATTN 

D.  Bernstein 

Kaman  Sciences  Corp 

ATTN:  Library  Science  Appl ications,  Inc 

ATTN:  W.  Layson 

Kaman  Tempo 


Southwest  Research  Institute 
ATTN:  A.  Wenzel 
ATTN:  W.  Baker 

SRI  International 

ATTN:  Y.  Gupta 
ATTN:  0.  Keough 
ATTN:  G.  Abrahamson 

Structural  Mechanics  Associates,  Inc 
ATTN:  R.  Kennedy 

Teiedyne  Brown  Engineering 
ATTN:  J.  Ford 
ATTN:  F.  Leopard 
ATTN:  D.  Ormond 

Terra  Tek,  Inc 


ATTN: 

S.  Green 

ATTN: 

J.  Schatz 

ATTN: 

Library 

ATTN: 

A.  Abou-Sayed 

Tetra  Tech  Systems  Co 
ATTN:  L.  Hwang 

TRW  Electronics  &  Defense  Sector 

ATTN:  Technical  Information  Center 
ATTN:  P.  Bhuta 
ATTN:  T.  Mazzola 
2  cy  ATTN:  N.  Lipner 

Universal  Analytics,  Inc 
ATTN:  £ .  Field 

Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  M.  Baron 

Weidlinger  Associates 

ATTN:  J.  Isenberg 

TRW  Electronics  &  Defense  Sector 
ATTN:  E.  Wong 
ATTN:  P.  Dai 


